Band filling and interband scattering effects in MgB2: C vs Al doping 
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We argue, based on band structure calculations and Eliashberg theory, that the observed decrease of Tc of Al 
and C doped MgB2 samples can be understood mainly in terms of a band filling effect due to the electron doping 
by Al and C. A simple scaling of the electron-phonon coupling constant A by the variation of the density of states 
as function of electron doping is sufficient to capture the experimentally observed behavior. Further, we also 
explain the long standing open question of the experimental observation of a nearly constant tt gap as function 
of doping by a compensation of the effect of band filling and interband scattering. Both effects together generate 
a nearly constant tt gap and shift the merging point of both gaps to higher doping concentrations, resolving the 
discrepancy between experiment and theoretical predictions based on interband scattering only. 



The high critical temperature of 40 K in the simple binary 
compound MgB2 (Ref. 1) was an unexpected present of na- 
ture to the scientific community. Now, after a few years of in- 
tense experimental and theoretical research the main features 
of superconductivity in this material seem well understood as 
due to a phonon mediated mechanism with different coupling 
strengths to the electronic a- and tt -bands , which 

leads to the appearance of two distinct superconducting gaps. 

Historically, two-band superconductivity is an old topic 
which has been proposed already shortly after the formulation 
of the BCS theory. Suhl, Matthias and Walker |6] suggested a 
model for superconductivity in transition metals considering 
overlapping s- and rf-bands. At the same time, Moskalenko 
formulated an extension of BCS theory for multiple bands ^^l- 
In the early 1960 there have been experimental claims for the 
observation of two-band superconductivity in some transition 
metals like e.g. V, Nb and Ta ^SJ and later, in the 1980, in 
oxygen depleted SrTiOs 01 ■ 

However until now, MgB2 appears to be first system for 
which multi-band superconductivity has independently been 
evidenced by several experimental techniques: heat capac- 
ity, tunneling spectroscopy, Raman spectroscopy, penetration 
depth measurements, ARPRES and the analysis of the critical 
fields |2J . The appearance of multiple gaps had been predicted 
theoretically |4] based on the electronic structure of MgB2 
The Fermi surface consists of four sheets: Two 
cylindrical sheets corresponding to quasi two-dimensional a- 
bands and two tubular networks derived from the more three 
dimensional tt -bands |5]. The phonons, in particular the op- 
tical bond-stretching phonon branch along F-A 1 11], couple 
about three times stronger to the holes at the top of the cr-band 
as compared to the 7r-band REHEEHESl- Using linear 
response theory it is possible to calculate from first principles 
the electron-phonon coupling (Eliashberg function) which is 
needed as input for Eliashberg theory. The solution of the 
Eliashberg equations allows for the calculation of the super- 
conducting gaps or thermodynamical properties like specific 
heat in good agreement with the experiments 



As for any anisotropic order parameter, scattering by non- 
magnetic impurities should have a pair breaking effect, just 
as magnetic impurities have in conventional superconductors. 
Interband impurity scattering should lead to a decrease of 
and if strong enough to a single (averaged) order parameter 
The interband impurity scattering between the a- 
and TT-bands is exceptionally small 1 18], due to the particular 
electronic structure of MgB2 , so that in the superconducting 
state, the two gaps in the a- and the 7r-bands are preserved 
even in 'dirty' samples with a considerably reduced and a 
broad range of normal state resistivities. 

The decrease of Tc has been experimentally demonstrated 
by a series of substitution experiments in which Mg has been 
replaced bvAl B hy C \^ 2^ 2L 21. 22^ 2A. 2^ 2^ 
|23,|28l|23,|3fl|3l|l- Similarly, irradiation with neutrons leads 
also to a decrease of Tc 13211 . Figure shows a compilation 
of experimental data for the critical temperature Tc versus Al 
and C doping concentration. 

For the two superconducting gaps it has been observed that 
the a gap decreases with decreasing Tc and approaches the in- 
termediate coupling value of 2 A Ik-Q Tc at Tc ^ 25 K. In most 
experimental reports the tt gap is found to be independent on 
the Tc of the sample and to remain close to the value of ^ 2 
meV seen for undoped samples |42]. 

There have been recent reports by Gonnelli and cowork- 
ers which demonstrate a different behavior of the 
superconducting gaps depending on the type of dopant. For 
C doped single crystals with composition Mg(Bi_j,Cy)2 
(j/<0.132) their point contact spectroscopy measurements 
show a merging of the a and tt gaps for the first time 1 33, 43]. 
On the other hand, Gonnelli et al. also find that the behavior of 
Al doping in single crystals with composition Mgi_2;Al2;B2 
{x < 0.21) is quite different. Even samples with very low Tc 
of about 20 K still exhibit distinct gaps, at critical tempera- 
tures for which theoretical calculations based on Eliashberg 
theory including interband scattering always predict a single 
order parameter only |44]. Therefore, at present, there is dis- 
agreement between experiment and theory. 
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FIG. 1: Critical temperature Tc as function of Al (filled sym- 
bols) and C (open symbols)versus doping concentration defined 
as Mgi_,Al,(Bi_^Cy)2. Al-doping: (m)M, (A)|28], (.)|26], 
(♦)|341, (►)1351 C-doping: (*)|361,^ (c>)|^ (^)| 30], (a)\37l 
(0)|38], (x)(3il|Q, (A)|40], (+)|29], (<l)(41]. The lines present 
estimates based on Eliashberg theory with different levels of approx- 
imation: dotted line with DOS from rigid band model, dashed line 
with DOS from VGA and solid line with DOS and phonon renormal- 
ization from VGA. The horizontal dotted line is the lower limit for 
interband scattering only. 



In the following we will argue that one essential ingredient 
to understand the behavior of Tc is the effect of band filling 
of holes in the cr-band due to electron doping. To understand 
the different behavior of the two gaps in Al and C doped sam- 
ples one additionally needs to consider interband scattering. 
While band filling will decrease the superconducting gaps, in- 
terband scattering will decrease the value of the larger gap and 
increase the smaller one. These two effects may compensate 
for the smaller tt gap and enable us to explain the observed 
nearly constant value of the small gap. 

First, we will focus on the doping dependence of the critical 
temperature. Figure [J summarizes experimental results from 
different groups. Tc as function of Al and C doping shows 
very similar behavior if the C doping is scaled by a factor 
of two as compared to the Al doping. This follows naturally 
from the definition of the C doping concentration per boron 
atom, as expressed in Mgi_a;AIa;(Bi_j,Cy)2, with x (y) for 
the amount of Al (C) doping. The importance of the band fill- 
ing is already indicated by the horizontal dotted line ^^25 K. 
This value would be the upper limit of Tc due to the pair break- 
ing effect of interband scattering only. If interband scattering 
would be the only relevant mechanism no sample should show 
a Tc lower than indicated by the horizontal line. This is clearly 
not the case as shown in Fig. [2 

MgB2 has a total of 0.26 holes: 0.15 holes in both cr bands 
and the remaining 0.11 holes in the hole 7r-band. Al and C 



substitution will both dope electrons and therefore reduce the 
number of holes. In a rigid band model the electron doping 
would be defined with respect to the total number of holes in 
MgB2 and simply corresponds to a shift of the Fermi level. 
For small doping the cr-band DOS is practically constant as 
expected from the quasi two-dimensional character of the a- 
bands. After adding 0.15 electrons the cr bands become nearly 
filled and the DOS starts to decrease rapidly. The coupling 
of the cr-holes to the optical bond-stretching E2g phonons 
drives the superconductivity in this material and determines 
Tc- Therefore, we just scale the band splitted electron-phonon 
Eliashberg functions a^jF [14] and the /i* -matrix 1451] with 
the the change of the a- or 7r-band DOS as function of dop- 
ing. We use the Eliashberg functions for pure MgB2 calcu- 
lated from first-principles linear response theory II Iv . which 
have been used successfully to describe the specific heat (3> 
tunneling |44|] and penetration depth 146,i] . 

The dotted line shown in Fig. [l] corresponds to the rigid 
band scaling. The decrease in Tc for small doping concentra- 
tions is well reproduced and originates from the small dis- 
persion of the cr bands along the F-A line. The cr-band Fermi 
surfaces are not perfect cylinders but are slightly warped (see 
Fig. 3, Ref. ]5]). For larger doping concentration, Tc obtained 
from this simple model decreases faster than observed in ex- 
periment. This is not surprising because we used the unper- 
turbed band structure of pure MgB2 not taking into account 
neither alterations of the bands due to doping nor the change 
of the phonon frequencies. 

To correct for this failure we further calculate the change 
of the DOS using the virtual crystal approximation (VCA). In 
order to simulate the doping of Al, we replace the Mg atom 
with a virtual atom with charge Z = xZai + (1 — x)ZMg and 
recalculate the electronic band structure self-consistently us- 
ing the full potential LMTO method ['TT]. In agreement with 
|j4SJ, we find a slower decrease of the cr-band DOS. Using 



the DOS from the VCA to scale a?, F we solve the Eliash- 
berg equations and obtain a slower decrease of Tc (dashed line 
in Fig.Q in better agreement with the experimental observa- 
tions. Recent supercell calculations indicate an even slower 
cr-band filling] 49] compared to the VCA. 

An additional effect of doping will be the hardening of the 
E2g phonon branch which will decrease the electron-phonon 
coupling A ^ l/uj"^ ] 10]. In order to take this effect into ac- 
count we also calculated the i?2g-r-point frequency in the 
VCA using linear response methods ITtI] . The final result 
from scaling af^Fhy the DOS and the E2g phonon frequency 
is shown by the solid line in Fig.^ The agreement with ex- 
periment improved significantly. 

Band filling with the corresponding changes in the DOS 
seems to be sufficient to understand the behavior of Tc as 
function of doping. However, this is not sufficient to under- 
stand the evolution of the superconducting gaps, because there 
should be no difference in the behavior between Al and C dop- 
ing because both are electron dopants. We now argue that the 
additional ingredient to understand this behavior is interband 
scattering. 
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FIG. 2: Upper panel Al doping: Superconducting cr-gap (upper curve 
and filled symbols) and vr-gap (lower curve and open symbols) as 
function of critical temperature Tc obtained from the solution of the 
Eliashberg equations with scaled afjF without interband scattering 
(solid lines) compared to experimental results (□) |33] and (o) |28]. 
Lower panel C doping: The solid lines show the solution of the 
Eliashberg equations with interband scattering rate 2000 ■ y crcT^ 
compared to the experimental results (□) IsslEsll . The limiting cases 
of interband scattering only (dotted lines) or scaling of a\jF (dash- 
dotted lines) are also shown. The dashed straight line indicates the 
BCS gap relation for A/Tc 



In the upper panel of Fig.|2we plot the experimental a- and 
TT-gaps for Al doped crystals as obtained by Gonnelli et al. 
13311 and by Putti et al. ilm as function of the critical temper- 
ature of the samples. Together with the experimental data we 
display the results from the solution of the two-band Eliash- 
berg equations without interband scattering but the Eliashberg 
functions scaled by the change of DOS and phonon frequency 
as described above. 

The agreement with experiment seems to be reasonable. 
The results by Putti et al. |28] show a merging of the gaps 
for Al-doping, however both gaps have a ratio of A/ZcsTc 
somewhat lower than the canonical BCS ratio as indicated by 
a dashed line in Fig.|2] which casts doubt on this data point. 
However, it is still premature to give a definite answer because 
more experimental data for single crystals and high doping 
concentrations will be required for a complete picture. Based 
on the available data which show no clear merging of the two 



gaps, we conclude that the interband scattering in Al doped 
samples is small, even at high doping concentrations. 

The experimental results for C doped single crystals indi- 
cate for the first time a merging of the two superconducting 
gaps, which is a clear manifestation of interband scattering. 
In the lower panel of Fig.|2lwe show the experimental results 
f33il together with our Eliashberg theory calculations. In dif- 
ference to the previous case we include the interband scatter- 
ing in our calculations, which will also cause an additional re- 
duction of Tc- Using a simple linear relation of the interband 
scattering rate to the doping concentration {1{^xq^ = 2000 • y 
cm^^) we find e.g. for 10-15% C concentration an additional 
lowering of Tc of about 6 K. 

The two limiting cases (interband scattering only, scaling 
of ajjF) are also shown in the lower panel in Fig.|2l The de- 
crease of the DOS causes a decrease of both gaps, as can be 
seen from the dash-dotted lines. In contrast the interband scat- 
tering will decrease the a gap and increase the tt gap. Both 
effects can compensate each other resulting in the solid line, 

which includes the effects of the scaled a? F and interband 

'■J 

scattering. This may explain the experimental observation of 
a nearly constant 7r-gap as function of doping, which has been 
a long standing open question. 

The difference in the magnitude of the interband scattering 
for Al and C doping can be easily rationahzed. The a bond- 
orbitals are located in the boron plane and there is not much 
weight of the cr-bands in the Mg plane. The 7r-orbitals are also 
centered at the boron plane, but extend further out towards the 
Mg plane. For that reason impurities in the boron plane are 
more effective interband scatterers 1 18, 50]. Therefore inter- 
band scattering due to doped C atoms replacing B atoms is 
much more likely than for Al doping. 

In summary, we have shown that the variation in Tc of Al 
and C doped samples of MgB2 can be understood mainly as 
due to a simple effect of band filling. Al and C are both elec- 
tron dopants which reduce the number of holes at the top of 
the (7 bands together with a reduction of the electronic DOS. 
Further, we suggest that the nearly constant tt gap as function 
of doping can be understood due to a compensation of band 
filling and interband scattering. The compensation of these 
effects shifts the merging point of both gaps to higher dop- 
ing concentrations and lower Tc, resolving the discrepancy 
between experiment and theoretical predictions based on in- 
terband scattering only. 
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